Optically guided beam splitter for propagating matter waves 
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We study experimentally and theoretically a beam splitter setup for guided atomic matter waves. 
The matter wave is a guided atom laser that can be tuned from quasi-monomode to a regime where 
many transverse modes are populated, and propagates in a horizontal dipole beam until it crosses 
another horizontal beam at 45°. We show that depending on the parameters of this X configuration, 
the atoms can all end up in one of the two beams (the system behaves as a perfect guide switch) , or 
be split between the four available channels (the system behaves as a beam splitter). The splitting 
regime results from a chaotic scattering dynamics. The existence of these different regimes turns 
out to be robust against small variations of the parameters of the system. From numerical studies, 
we also propose a scheme that provides a robust and controlled beam splitter in two channels only. 

PACS numbers: 37.27+k,37.10.Gh,05.45.-a,67.85.-d 



Matter- wave interferometry has the potential of being 
several orders of magnitude more sensitive than optical 
interferometry. Many demonstrations have been made of 
its use as high-precision gravimeters or gyrometers [1]. 
Its development requires the coherent manipulation of 
matter waves with atom optical elements. In this respect, 
the achievement of beam splitters for atoms moving in 
free space was essential [2-8]. 

To optimize the sensitivity of an atomic interferome- 
ter, one needs large angle beam splitters with atom wave 
packets into superposition of two narrow momentum dis- 
tribution with large momentum separation. Narrow dis- 
tribution yields good fringe contrast [9]. A few strate- 
gies are currently explored or envisioned to increase the 
enclosed area: the enhancement of the interaction time 
based on large momentum beam splitter ^7|, i8|, ilO|] , atom 
interferometry experiments placed in a reduced gravity 
environment [ijj, ll2| , or the use of slow atoms in a guided 
environment. 

In this Letter, we shed light on the physics of matter 
wave splitting in a crossing guide configuration. Con- 
fined geometry for atom interferometry are promising in 
terms of compactness and portability. The wide variety 
of techniques available to design potentials and guides 
for the external degrees of freedom has been used to in- 
vestigate the operation mode of different type of beam 
splitters. One should distinguish those involving a trap 
[IQI , [l 3| -[ l 8l and those based on combinations of waveg- 
uides [19'- 23]. So far, the former have been investigated 
with Bose-Einstein condensate (BEC) while the latter 
were studied in the thermal regime. In this work we 
report on the study of an all-optical beam splitter struc- 
ture for propagating matter waves both experimentally 
and numerically and in both regimes (i.e. various ini- 
tial transverse modes of the guide from the multimode 
to the monomode limit [24]). This problem amounts to a 



quantum scattering problem in confined environment in 
which chaotic behavior can emerge J25l-[27j. 

The beam splitter is obtained by combining two guides 
which provide a potential structure where four paths are 
available for the matter wave. In our experimental setup, 
we use two dipole beams crossing at 45° in an horizontal 
plane (see Fig.[T]). This X configuration has four channels 
[28]. Optical confinement has many advantages: it can be 
applied to atoms with no permanent magnetic moment 
(e.g., ytterbium and alkaline-earth), which is particularly 
interesting for metrology, an extra magnetic field can be 
used to tune the interactions (Feshbach resonance) and 
the focus position can be easily moved [29|, [30|]. Fur- 
thermore, the crossing region has a confinement strength 
larger than that of the guides. As a result the setup is 
robust against low frequency noise. This is to be con- 
trasted with other Y geometries where the connection 
zone between the guides implies an important lowering 
of the potential strength [19~21i , i31i] . 

We produce a rubidium-87 BEC in a crossed dipole 
trap (located at A in Fig. [1]) using two focussed laser 
beams (wavelength of 1070 nm, waist of wi c^ 35 jam) 
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FIG. 1: (color online) Schematic of the matter wave guided 
beam splitter setup based on two crossing dipole beams in a 
horizontal plane. Atoms are out coupled from a BEC located 
at point A and propagate in the beam Li towards point B 
where beam L2 crosses beam Li. 



[245. Beams Li and L2 originate from the same laser 
but have a frequency difference of 80 MHz to wash 
out possible interferences. Atoms are prepared in the 
\F = l^rriF = 0) internal state through a spin distilla- 
tion process implemented during the evaporative cooling 
stage [32| . The guided atom lasers are realized by outcou- 
pling atoms from the BEC in the horizontal dipole beam 
of the trap either by applying a time-dependent magnetic 
field gradient or by reducing the intensity of the non hor- 
izontal beam of the crossing dipole trap [2J, |32| - |36| . 

Although atom-atom interactions may not always al- 
ter dramatically the coherence [37], theoretical studies of 
cold atoms interferometers based on guiding potentials 
show that they can affect dramatically the contrast of the 
expected interference signal [38|, |39| . To reduce the role 
played by interactions we perform the outcoupling very 
slowly. The low atomic density of the propagating beam 
makes this system nearly interaction- free [24, 32]. The 
propagating matter wave is charaterized by its atom flux 
(~ few 10^ atoms/s), its mean velocity (13 ± 2 mm.s~^) 
and the mean excitation number (n) associated to the 
transverse modes. 

Experiments are carried out at a fixed power Pi = 
140 mW of the guide beam Li which corresponds to a 
measured transverse frequency of about 300 Hz and a 
depth of 10 /iK. The crossing beam L2 has a waist W2 = 
70±5 /im and an adjustable power P2. Atoms propagate 
for a time long enough to obtain clear information on 
the asymptotic output channels. In practice, the crossing 
between Li and L2 beams is at 700 jam. downstream from 
the BEC. The beam propagates during 200 ms, and we 
take an absorption picture after a few ms time-of- flight. 

In a first set of experiments, we prepare the prop- 
agating guided atom laser in the ground state of the 
transverse confinement ((n) !::^ 0). In this transverse 
mo no mode regime, we observe three different regimes 
(I), (II) and (HI) depending on P2/P1 (see Fig. [2|). Im- 
ages (I) belong to the weak defect regime for which the 
beam is transversally excited after passing through the 
interaction region. The coherent excitation of the trans- 
verse modes is induced by the weak coupling between 
the longitudinal and transverse degrees of freedom in the 
crossing guide region. As illlustrated in Fig. [2l the larger 
P2/ Pi^ the stronger the coupling and the larger the num- 
ber of transverse modes populated downstream from the 
crossing region. Image (II) corresponds to the splitting 
regime where atoms can explore the four channels pro- 
vided by our X shape geometry. It corresponds to a much 
stronger coupling between the different degrees of free- 
dom and the correspon ding classical dynamics becomes 
chaotic (see below) [23, |40| . Images (HI) illustrate the 
switch regime where atoms go into one arm of the beam 
L2 with a 100 % efficiency [20|, [21|. The center of mass 
motion is excited as a result of momentum transfer of 
the longitudinal incident momentum on the transverse 
degrees of freedom of the output channel guide. We plot 
a boundary between (I) and (II) in Fig. [2] (middle) when 
less than 80 % of the atoms continue their propagation 



in the original guide Li, and between (II) and (HI) when 
more than 80 % of the atoms are deviated to the guide 
L2. The error bars show the typical size of the transition 
regions. 

We also have investigated experimentally this dynam- 
ics for various values of the initial mean quantum num- 
ber (n) of the transverse states by adjusting the frac- 
tion of condensed atoms from which the guided atom 
laser is produced [2J]. For (n) = 10 more than hundreds 
of transverse levels are populated since one has to take 
into account the degeneracy, while one transverse state 
is populated for (n) =0. The results are summarized in 
Fig. [2] (middle). We observe that the three regimes are 
still present. However, the transverse oscillations can be 
observed only for low incoming mean transverse excita- 
tion (compare images (Ilia) and (Illb) taken for different 
values of (n)). In previous experiments where the split- 
ting of a beam of cold atoms was carried out, it was im- 
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FIG. 2: (color online) Propagation of an atom laser in a X 
shape guide structure. Top: Three regimes are experimentally 
observed depending on the power ratio associated with each 
guide: (I) weak defect regime where the excitation of trans- 
verse modes downstream from the crossing increases with the 
power ratio P2/ Pi^ (II) splitting regime where atoms explored 
all the available channels, and (Illa/b) switch regime where all 
atoms are directed into guide L2. The mean transverse exci- 
tation numbers of the incoming beam are (n) :^ 0.3 (Illa) and 
(n) :^ 10 (Illb). Middle: Experimental diagram that summa- 
rizes the observed regime depending on both the power ratio 
P2/P1 and the mean transverse populations. The boundaries 
of the three regimes turn out to be robust against the popula- 
tion of the transverse modes. Bottom: Numerical simulations 
confirming this robustness; red solid line: quantum simula- 
tions, blue dotted line: classical simulations. Wave packets 
are started 150 /im before the crossing, with initial longitudi- 
nal velocity 10 mm.s~^, for w\ = 45 //m, W2 — 41 /xm. 



possible to observe this transverse excitation since many 
transverse states were populated and these oscillations 
were thus washed out [l9|, [20|, S, 123 • The boundaries 
between the different regimes appear to be very robust 
against (n) and a slight misalignement of the crossing 
beam. This is to be contrasted with the interaction of a 
guided atom laser with an out-of-center vertical poten- 
tial [27| . The boundaries remain nearly unchanged when 
the crossing is misaligned by a tenth of the waist of the 
beams. For each set of data the alignment was checked 
before and after the data acquisition. We estimate that 
the maximum relative change of the beam positions is 
below 2 jam. 

To interpret the splitting regime, we have computed 
the dynamics of 2D quantum wavepackets using the split- 
Fourier algorithm [41^ and developed a direct simulation 
Monte Carlo (DSMC) method where atoms are evolved 
according to classical mechanics (see Fig. [2] (bottom)). 
The three regimes observed experimentally can be re- 
produced by numerical simulations and turn out to be 
generic. Each point in Fig. [2] (bottom) corresponds to 
the results obtained for an initial state that coincides 
with the n}^ transverse eigenstate. We observe that the 
boundaries deduced from classical and quantum simula- 
tions are very close. The fact that the boundaries are 
not exactly the same in classical and quantum dynam- 
ics is somehow expected since the equations fulfilled by 
the center of mass of the packet are different in classical 
and quantum mechanics. Indeed, the external potential 
is not anymore harmonic in the crossing guide region. 
The error bars on the experimental data are such that 
one cannot distinguish the classical or quantum nature 
of the motion at low n. Numerics confirm the robustness 
of the existence of the boundaries. However, their precise 
positions can depend on several parameters such as small 
misalignment of the crossing beams, a residual longitu- 
dinal acceleration and/or the relative waist size of the 
two beams. Many of these parameters cannot be mea- 
sured accurately. Consequently and to avoid the use of 
too many free parameters, we do not fit the experimental 
data with the numerics but use the numerics to identify 
the underlying physics and confirm its generic character 
and its robustness. 

Figure [3]provides a typical experimental example of the 
dynamics of the splitting regime. The guided atom laser 
that enters the crossing region remains partially trapped 
(larger density) and then the contamination of all avail- 
able branches is observed. We have plotted in Fig. 3] 
the output channels obtained from a simulation relying 
on classical mechanics (using different color for the differ- 
ent output channels) for different initial conditions of the 
transverse phase space. We clearly identify a chaotic re- 
gion with a fractal structure in which a very small change 
in the initial condition changes the output channel. This 
is the signature of a chaotic dynamics [25]. We have also 
checked that Lyapunov exponents become positive in this 
parameter region. To further characterize this zone, we 
have verified that the set of boundaries between the four 



accessible output channels as a function of the initial con- 
ditions displays a fractal structure. For a wide variety of 
fixed transverse velocity Vy in regime (II), we find a frac- 
tal dimension of ^ c^ 0.85 revealing a strong fractality 
of the basin boundaries, usually associated with chaotic 
scattering. Furthermore, numerics seem to indicate that 
the fractal Wada property [25|, |42|, |43| (every point on the 
boundary of a basin is on the boundary of all basins) can 
be fulfilled in this system. The distribution of trapping 
times in the vicinity of the crossing guide region is found 
to exhibit an exponential decay. 

Interestingly the population in a given channel can be 
reduced by adding a small extra gradient ; this happens in 
the experiment if the beams are not perfectly horizontal. 
In our nearly horizontal configuration, a gradient that 
corresponds to a hundredth of the gravity strength is 
sufficient to dissymmetrize dramatically the population 
in the output channels. In this way one can populate only 
two channels and realize a two-channel beam splitter. 
This effect has been checked numerically and explains 
the variation of the populations in the different channels 
seen in the experimental results (compare image (II) of 
Fig. [2] and images of Fig. [3]). A beam splitter experiment 
for a thermal cloud accelerated by gravity in a quasi- 
vertical configuration has been discussed in [22|, Q . Note 
that the possibility of realizing a coherent beam splitter 
using the transient dynamics of a chaotic system has been 
addressed in another context with amplitude modulated 
optical lattices [45]. 

The X configuration that we have investigated thus 
exhibits two regular regions separated by a chaotic re- 
gion. Actually, the numerical simulations reveal the ex- 
istence of other splitting region at even higher power ra- 
tio P2/P1 that were not investigated in our experiment. 
These splitting regimes are undesirable for interferome- 
try since the relative phases of the matter wave in the 
different output channels can be affected in an uncon- 
trollable manner from run to run. A possible solution as 
investigated numerically in [46] consists in adding a ded- 
icated extra potential to decrease the coupling between 
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FIG. 3: (color online) Example of experimental splitting dy- 
namics illustrating regime (II) (see Fig. [2|). Time evolution 
is directly visible in the succession of absorption images. The 
images show a larger density at the crossing point which is at- 
tributed to the complex dynamics that traps the atoms there 
for some time. Atoms that escape from this region populate 
all available channels. 



the longitudinal and transverse degrees of freedom. 

In the following, we propose an alternative scheme us- 
ing a transverse momentum kick transferred to atoms 
before they enter the crossing region. The power ratio is 
chosen to be in regime (III) (the switch guide regime). 
Different experimental techniques can be used to achieve 
such a transfer: Bragg or Raman transitions, Bloch oscil- 
lations [47|. In our simulation, we add an instantaneous 
transverse momentum upstream to account for the mo- 
mentum kick. Figure [5] summarizes our numerical results. 
Each asymptotic output channel obtained here from a 
classical simulation is represented by a different color as 
a function of the velocity kick value and the "phase" i.e. 
the position before the crossing at which the momen- 
tum transfer is performed. Using the prediction of the 
simulations relying on classical physics as guidelines for 
quantum simulations, we conclude on the feasibility of 
re-directing atoms in guide Li with an efficiency larger 
than 96 %. This technique enables to switch from guide 
1/2 to Li without changing in time the guide characteris- 
tics [20|, l21| . This also means that if the incoming wave 
packet is prepared in a linear superposition with equal 
weight of the states with and without the proper trans- 
verse momentum kick, the wave function can be split 
into two coherent wave packets that propagate in the two 
guides Li and L2. Interestingly enough, the domain of 
useful parameters in momentum kick and phase to real- 
ize such a guided matter wave beam splitter is relatively 
large. This shall ensure the robustness of the method. 
The splitted beams are for each guide in an excited state 
that is responsible for their transverse oscillation (see im- 
age (III) of Fig. [2j). An out-of-guide axis potential with 
an appropriate strength can be used to remove the re- 
maining transverse dipole oscillation [27]. 

In conclusion, this Letter reports on the realization of 
all-optical guided beam splitter devices. Such tools are 
important for the development of guided atomic optics 
which has a huge potential for applications, in particular 
for inertial sensors. In addition, it was shown recently 
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FIG. 5: (color online). Left: Output channels for a classical 
trajectory started at the bottom of the guide, between 100 
/im and 200 /xm before the crossing, with initial longitudinal 
velocity 10 mm.s~^, and an initial transverse velocity between 
—25 mm.s~^ and 25 mm.s~^, with wi — 45 /xm, W2 — 41 
/im, P2 — 0.365 W= 2.04 Pi. Right: center of mass motion 
of numerically simulated quantum wave packets. The one 
with initial parameters of point C (150 /xm before defect, no 
transverse velocity) is deviated in the second guide at 99.9 % 
(green curve). In contrast, the one with initial parameters of 
point D (154.75 /xm before defect, transverse velocity of 15.5 
mm.s~^) ends up at 96.2 % in the initial beam (blue curve). 



that there exists a direct mapping between waveguide 
theory and spin chain transport [48|, opening possible 
applications of our scheme in solid state physics. 



FIG. 4: (color online). Left: colors associated with the dif- 
ferent output channels. The arrow indicates the propagation 
direction of the incoming atoms. The overlap region between 
the guide is represented as a dark blue disk. Right: Output 
channels associated to different initial phase space conditions 
of the transverse degree of freedom in regime (II) (see Fig.[2|). 
Each classical trajectory started at the bottom of the guide, 
150 /im before the crossing, with initial velocity 10 mm.s~^, 
for wi = 45 /im, W2 = 41 /im, P2 = 0.143 W= 0.8 Pi. 
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